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A PHENOMENOLOGICAL RELATION BETWEEN STRESS, 
STRAIN RATE, AND TEMPERATURE FOR METALS 
AT ELEVATED TEMPERATURES 
By Elbridge Z. Stowell 

SUMMARY 


A phenomenological relation between stress, strain rate, and temper- 
ature is suggested to account for the behavior of poly crystalline metals 
above the equicohesive temperature. The properties of the metal 
included in the relation are elasticity, linear thermal expansion, and 
viscosity. The relation may be integrated under various conditions to 
provide information on creep rates, creep rupture, stress-strain curves, 
and rapid-heating curves. It Is shown that for one material - 7075 -T6 
aluminum-alloy sheet - the information yielded by the relation for these 
four applications agrees reasonably well with test data. 


INTRODUCTION 


The behavior of metals at elevated temperatures is constantly 
becoming a more important problem. Not only must information be at 
hand concerning creep at these temperatures, but Information must also 
be available concerning when rupture might be expected if the creep 
were allowed to continue, what the stress-strain curve is likely to be 
at some particular temperature and strain rate, and how the material 
will behave when rapidly heated to a high temperature. It would be a 
useful achievement if some fundamental- law characteristic of the material 
could be found, such that these and other apparently distinct problems 
could be shown to be simply different aspects of the application of one 
law. 


In studies of the creep of metals, it has become apparent that the 
creep behavior at elevated temperatures is different from that at low 
temperature s . For relatively pure metals, the line of demarcation for 
temperature occurs at about one-half the melting point and may be tenta- 
tively identified with the equicohesive temperature. In the high- 
temperature region, certain laws have been found and explicitly stated 
(ref. l). These laws appear to give some promise that the behavior of 
metals in this region may be predicted for effects other than creep. 
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Creep is characterized by an initial transient or primary region 
which is followed by a steady or secondary region in which viscous flow 
occurs and finally by an accelerating or tertiary region which leads to 
rupture . A theory based on rate -process considerations is already 
available and accounts in a qualitative manner for the main effects 
shown by the steady creep in the high-temperature region (ref. 2). For 
stable metals, this theory requires a knowledge of only three constants 
to specify the creep behavior: one constant is roughly predictable from 

the melting point, and the other two must be determined from tests. 

In this report, a phenomenological relation between stress, strain 
rate, and temperature is suggested for this region, and this relation 
utilizes the existing rate-process theory for steady creep. The rela- 
tion is not an equation of state of the type proposed by Hollomon 
(ref. 3) but is an equation of "rate of change of state." Since dif- 
ferent integrals of this relation exist under different conditions, 
there is no single "equation of state" in this case. These integrals 
are applied to creep, creep rupture, stress-strain curves, and rapid- 
heating curves, and comparisons of these numerical results are made with 
experimental results for 7075 -T 6 aluminum- alloy sheet. 


SYMBOLS 


€ strain 

t time, hr 

e strain rate, per hr 

€ Q strain rate used in obtaining stress-strain curves, per hr 

a stress, ksi 

a Q constant, ksi 

CT lim limiting stress, ksi 

E Young ' s modulus , ksi 

a linear expansion coefficient, per °K 

T temperature, °K unless otherwise indicated 

T rate of temperature change, °K per hr unless otherwise 

indicated 
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AT = T 0 t 

AFT activation energy, cal per mole 

R gas constant, taken as 2 cal per mole per °K 

b constant, per hr per °K 


. & 
= Jo. ~RT 

2sT 


H 


= VI + 


2p 



Subscript: 

r conditions at rupture 


STATEMENT OF PHENOMENOLOGICAL RELATION 


A metal undoubtedly retains a certain elasticity at elevated tem- 
peratures, even though the elastic range is reduced from that prevailing 
at lower temperatures . The elasticity resides in the crystalline grains 
and is a measure of the resistance of the crystal lattice to distortion 
by stress. In addition to the elasticity, a polycrystalline metal also 
possesses a quasi-viscosity in the high-temperature region; in this 
region the application of a constant stress will result in a constant 
strain rate. Also, thermal expansion (or contraction) due to temper- 
ature change will create an apparent strain which can contribute to the 
total strain. 

The suggested phenomenological relation which includes all these 
properties and which takes into account possible temperature changes 
may be stated ' as follows : 


_£E 

M = k = A- (£) + a SC + 2sTe 51 sinh -f- (l) 

dt dt \E/ dt <*o 

In this relation the first term on the right supplies the contribution 
to the strain rate due to elasticity; the middle term, the contribution 
due to 'possible thermal, expansion; and the final term, the contribution 



4 


NACA TN 4000 


due to viscosity. 'A relation somewhat similar to this "but without the 
middle term has been proposed for plastics. (See ref. 4 .) 


The form of the viscous term is taken from reference 2 and is now 
well established. The exponential term describes the principal effect 
of temperature. The hyperbolic sine term describes the effect of stress. 


Kauzmann in reference 2 gave this term as 


sinh 



where 


m 


is 


a constant, but it has been demonstrated repeatedly by tests that the 
argument of the term is independent of temperature and is of the form 


— , where a is a constant. The remainder of the expression, which 


consists of a number of individual constants multiplied by the temper- 
ature T as given in reference 2, is replaced here by 2sTj the factor 2 
is retained for convenience. 


The quantities E and a are values which are usually known 
as functions of the temperature . The quantities s, Al/R, and a Q 
are found from measurements of the steady creep of the metal under 
constant stress and temperature. A range of temperature in the high- 
temperature region and also a range of stress are required to establish 
AH/ R and cr 0 , respectively. 


APPLICATIONS 


With different conditions imposed upon the metal, the suggested 
relation should generate in turn the laws governing such phenomena as 
steady creep, creep rupture, stress-strain characteristics, and behavior 
under rapid-heating conditions. These laws are derived from equation (l). 


Steady Creep 

If stress a and temperature T are constant, equation (l) 
reduces to 


■ 

e 


2sTe 


AH 

RT 


sinh 4L 
On 


( 2 ) 


or 


AI £_ 

RT Jo 
e 


€ = sTe 


( cr » q 0 ) 


(2a) 
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which indicates steady creep. Equation (2a) may also be written as 

-2- = ^ + log e 4 - log e sT (3) 

in which log g sT is substantially constant over a limited temperature 
range. Stress plotted against log g e should, therefore, yield straight 
lines with slope cr 0 which are displaced from, one another at different 
temperatures by amounts depending mainly on the value of Stress is 

sometimes plotted against the parameter ^ + logg e (ref. 5) to demon- 
strate that data taken at different temperatures will fall on a single 
curve. 


Creep Rupture 

At constant stress and temperature, the strain at a time t after 
the application of stress is 


e = It 

since the creep is steady under these conditions. 

Reference 6 discloses as an experimental fact that if t r is the 
time to rupture under conditions of creep, then 



» Constant = e r 


Thus, 



and equation (3) may be rewritten as 


g r _ AS 
CT o RT ~ 


log e t r + log 



e sT 


(3a) 


Since the last term is substantially constant, the rupture stress a r 

ATT 

plotted against the parameter ^ - log g t r should yield a single 
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straight line with slope a Q . Such a parameter has been employed in 
the correlation of creep-rupture tests (ref. 7)> 


Stress -Strain Curves 

At constant temperature, equation (l) becomes 

_ m 

e = i £2 + 2sTe" sinh -2- ; (4) 

E dt a Q v 

Obviously, if the stress is applied very rapidly, the effect of the 
final viscous term may be negligible compared with the 'term ^ ^2; ’ thus , 
the material may appear to be elastic. 

If an attempt is made to hold the strain rate constant, as is some- 
times done with conventional testing machines, by setting the strain 
rate equal to the testing rate e Q , a differential equation in a results, 
and the solution which is a family of theoretical stress-strain curves 
is given as follows: 


where 


« log. 


\|l + p 2 + (1 + p)tanh M 


\|l + p 2 + (1 - P)tanh M 


P = f°- e RT 
p 2sT 


M = \|1 + Ee 


2P a c 


e = e Q t 


(5) 
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Replacement of tanh M by exponentials gives 


— = log e 
ff o e 

and since ordinarily 
curves becomes 


(Vl + P 2 + 1 + &)e m + \j 

'l + p 2 - 1 - p 

Wl + p 2 + 1 - 3) e 24 + \ 

1 + p 2 - 1 + 0 


p » 1, the expression for the stress-strain 



(2p + ije 2 ** - 1 
e^ +2(3-1 


( 6 ) 


At the lower end of the curves where the strain is small, M is small, 
» i + at, and 


JL = log 2P + 2M + 
g Q & e 2p + 2M 


= log e 



2M6 \ 
M + Py 


~ agP- 

M + P 



which shows that the initial portion of the curves is elastic. At the 
upper end of the curves where the strains are large, M is large, and 
from equation (6), 


= log e (2p + 1) « log e 2P (7) 

where a-] is the limiting stress which can be developed at that tem- 
perature and strain rate. 

As a corollary, note that the following equation results when the 
value of 2p is substituted into equation (7): 


e Q = sTe 


jAH °llm 
RT e a o 


which is of the same form as equation (2a) for steady creep. Thus, the 
limiting stress ^ obtainable bears the same relation to the testing 

rate e Q that any creep stress g bears to its corresponding creep 
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rate e. This correlation has been noted experi m entally by Sherby and 
Dorn in reference 5 where the ultimate strength in tension is used for 

a lim‘ 


Rapid-Heating Curves 

The material composing the skin of a high-speed airplane or missile 
may be subjected to rapid increases in temperature as a result of aero- 
dynamic heating during flight. It is Imp ortant to know what to expect 
from the material under such conditions. 

In equation (1), let ~ = T_ in which T is a constant. Then, 

dt u ° 

for a constant stress a » a Q , 

J2_ AI 

e dt = 0 J d^ + aT 0 t + se CT ° J Te- RT dt (8) 

in which allowance should be made for the variation of E and a with 
temperature. Equation (8) may be written in terms of finite increments 
as follows: 

_ ^ 

Te RT At (9) 


= 3 
E 


+ a AT + se 


a 

au 


y 



where 


AT = T 0 t 

The initial term — gives the contribution to the total strain of the 
E 

application of stress, the second term gives the contribution due to 
thermal expansion, and the last term gives the contribution due to 
viscosity. 


COMPARISON OE THEORY AND EXPERIMENT 


A comparison of the laws generated by the assumed relation of 
equation (l) with experiments performed on one material, 7075 -T6 aluminum- 
alloy sheet, is given in this section. 
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The alloy is in the solution heat treated and aged condition. 
Further exposure to elevated temperature results in overaging. The 
material becomes unstable and suffers a marked reduction in strength 
with time when subsequently heated in the range from 3 00° to 600° F. In 
addition, different batches of the material may vary somewhat among 
themselves in properties. In view of this variability in the material, 
the constants in the viscous term were adjusted in each case to give a 
reasonable fit to the data. Table I shows the constants used in the 
calculations . 

The temperature which separates the high and the low regions is 
about 300° F for this alloy. Tests performed at temperatures below 
300° F should not be expected to show any agreement with calculations. 

In most cases , the results of such tests have been included to show the 
extent of the disagreement. 


Steady Creep 

Figure l(a) shows steady- creep data for 7075 -T 6 aluminum-all oy 
sheet. These data are obtained from reference 8 (solid curves ) and from 
preliminary unpublished NACA data (points). The dashed lines were 
computed from equation ( 3 ) by using the appropriate constants in 
table I. The dashed lines show a good fit to the data from reference 8 
at 300° F and 375° F; the disagreement at 211° F is plainly evident. 

The dashed lines computed for 450° F and 500° F go through the points 
obtained at the higher stresses in each case; whereas, the lower stresses 
give somewhat less creep than would be expected from equation ( 3 ). This 
peculiarity at low stresses has been observed previously by Dorn anri 
Shepard (ref. l). 

Figure 1(b) shows the same data plotted against the par ame ter 

— + log„ e. All the data fit reasonably well with the calculated curve 
RT e 

(eq. ( 3 )) with the exception of the points obtained at 211° F. 


Creep Rupture 

Figure 2 shows the creep-rupture data, taken from reference 8 
from preliminary unpublished NACA data, plotted against the parameter 
AFT 

— “ l°g e t r , where tj. is the rupture time at the stress cr r and the 

temperature T. Data show that a reasonable value for the constant e r 

Is 0.0077; the straight line computed from equation (3a), which uses this 
value and the values given in table I, goes through the lin e of points, 
with the exception as noted before of the very lowest stresses and of 
the points at 212° F. 
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Stress -Strain Curves 

Several investigators in the past have pointed cut relations 
between creep and stress-strain dataj for example, Hollomon (ref. 3) 
suggested a graphical procedure for passing from one to the other. In 
the present paper, the passage is accomplished analytically through equa- 
tion (6). 

Figure 3 shows stress-strain curves for the 7075 -T 6 aluminum- alloy 
sheet, as given in reference 9 , for three temperatures. The strain rate 
was 0.12 per hr. The theoretical curves obtained by computation from 
equation (6) (see table I) fit the data well at the higher temperatures 
but overshoot the 200° F curve by a considerable amount. 


Rapid-Heating Curves 

Figure 4 shows the rapid-heating curves for the 7075-16 al.um1.num 
alloy, t aken from reference 9 , for three stresses and for widely differ- 
ent temperature rates. The calculated curves represent a severe test 
of the assumed relation given in equation (l) not only because the 
cooperation of all four terms is required in equation (9) , which is 
used for computation, but also because the temperatures reached may be 
considerably higher and the strain rates considerably larger than those 
encountered in other types of tests. Near the upper ends of the curves, 
calculations were made at sufficiently close intervals of temperature, 
usually 5° to 1Q° K, to establish better the points. The shape and the 
separation of the curves with heating rate are given reasonably well by 
the theoretical curves except near 300° F, and then the predicted final 
rapid increases in strain occur at a somewhat higher temperature than 
that actually observed. The agreement at the lower temperature rate on 
the curve for 17-2 ksi Is not as good as the others. 


DISCUSSION 


It would appear from the tests reported here that the validity of 
the relation suggested in equation (l) Is borne out for the 7075-T6 alumi- 
num alloy. The relation has shown itself capable of accounting for the 
steady creep, the creep -rupture data, the stress -strain curves, and the 
rapid-heating curves for this alloy. 

A more satisfactory test of the validity of the relation would be 
with a more stable material such as pure aluminum. The alloy selected 
in this paper is not such a material but was used because of the exist- 
ence of the different types of data available. The elevated- temperature 
phenomena are reasonably well accounted for, however, in spite of this 
handicap. 
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The constants in the viscous term vary to different degrees. The 
activation energy £E of the three constants is most nearly constant 
and is structurally insensitive, being independent of stress, temperature, 
grain size, and minor alloy additions, and has been tentatively identi- 
fied with the activation energy for self -diffusion (ref. l). The magni- 
tude of ZsH is roughly proportional to the melting point of the metal. 

The quantity a Q} on the other hand, is structurally sensitive) for pure 

aluminum, it has a magnitude of about 0.2 ksi, whereas for the 7075 -T 6 
alloy it has a magnitude of k-.J ksi. For stable materials, it is 
independent of stress, temperature, and grain size. The quantity s 
was found to be a constant that is independent of stress and temperature 
and is subject to the same degree of uncertainty as the strain rate. 
Repeated steady-creep tests on the same material under identical con- 
ditions will yield creep rates within a range of perhaps 2:1. This 
uncertainty in the creep rate is reflected in a similar uncertainty in 
the quantity s. It is, therefore, not surprising to find a s imilar 
variation in s in the different applications. (See table I.) 

For materials which are unstable in the sense that the constants 
in the viscous term are mainly functions of the temperature, one might 
expect equation (l) to be still usable if such variations in the 
constants are allowed for in the computations. If the materials are 
unstable in the sense that the constants in the viscous term are mainly 
functions of the time, however, then equation (l) is probably not usable. 

Equation (1) is evidently valid for values of stress and temperature 
which are changing smoothly in one direction, even though these changes 
may be fairly rapid. For sudden or discontinuous changes in stress or 
temperature, transients are known to occur which will not be given by 
equation (l). Some modification of this relation, perhaps by an addi- 
tional term, will be required to account properly for such transients. 

Such a modification would have to account for transient creep as well 
as for other conditions; the nature of the modification is obscure at the 
present time. 

Certain parameters have appeared as a result of this relation and 
may prove useful in comparing data at different temperatures above the 
equicohesive temperature. For example, if steady-creep data are avail- 
able at different stresses and temperatures , stress plotted against the 

parameter ^ + log e k should yield a straight line with slope equal 
to a 0 . For the stress -strain curves, the parameter which determines 

. m 

the separation of the curves with temperature is 2^ = ~ e^ 1 , since 

S X 

the maximum height of the curves is cr 0 log e 2{3 . For rupture resulting 
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from creep, rupture stress plotted against ^ - logg t r should yield a 
single straight line with slope equal to cr 0 . If the rupture stress were 
plotted against the parameter — + log - log„ tp, the straight 

c S J. ^ 

line should pass through the origin. 

CONCLUDING REMARKS 


The application of the suggested phenomenological rate relation 
to the determination of steady creep rates, creep-rupture time, stress- 
strain curves, and rapid-heating curves for 7075-16 aluminum- alloy sheet 
has shown that such a relation can he used to cover a wide variety of 
applications. The possibility of predicting the behavior of materials 
under widely different load and temperature conditions is indicated for 
temperatures above the equi cohesive temperature for the material. Addi- 
tional verification of the validity of the relation would be desirable 
as suitable data over a sufficiently wide range of test conditions 
become available. 


Langley Aeronautical. Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 18, 1957- 
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TABLE I 

VAIDES OF THE CONSTANTS USED IN THE APPLICATIONS 



H 

■f=- 




10-3 


e, par Br 


(a) Creep rates at different temperatures • 


Figure i.- steady creep rates for 7075-Tb aJ.u mln um- a lj.oy sheet. 



Streesj a, k.si 
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d ) uorrexation or creep rates witn. parameter — +• log e e . 


Figure 1,- Concluded 
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Rupture stress. 



Figure 2.- Correlation of stress with rupture time for creep rupture of 

7075 -T 6 aluMnum-alloy sheet. 
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Figure 4 .- Rapid-beating curves for 7075-T6 aluminum-alloy sheet for 
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